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Abstract. - Mesoscopic aging systems are characterized by large intermittent noise fluc- 
tuations. In a record dynamics scenario [P. Sibani and J. Dall, Europhys. Lett. 64, 2003] 
these events, quakes, are treated as a Poisson process with average Qln(l + t/tw), where t is 
the observation time, is the age and a is a parameter. Assuming for simplicity that quakes 
constitute the only source of de-correlation, we present a model for the probability density func- 
tion (PDF) of the configuration autocorrelation function. Beside a, the model has the average 
quake size 1/q as a parameter. The model autocorrelation PDF has a Gumbel-like shape, which 
approaches a Gaussian for large t/t^n and becomes sharply peaked in the thermodynamic limit. 
Its average and variance, which are given analytically, depend on t/t^ as a power-law and a 
power-law with a logarithmic correction, respectively. Most predictions are in good agreement 
with data from the literature and with the simulations of the Edwards- Anderson spin glass 
carried out as a test. 

Introduction. - After a rapid quench of an external parameter, e.g. the temperature, 
many complex materials age, i.e. their properties slowly change with the waiting time, t^,, 
elapsed from the quench. Ever since the initial observations in polymers [1], evidence has 
accumulated that spin-glasses [2], type II superconductors [3], glasses [4], and soft condensed 
matter [5], among others, age in similar ways, e.g. : For observation times t <^ physical 
averages are nearly constant, and autocorrelations and their conjugate linear response func- 
tions are connected by an equilibrium- like fluctuation-dissipation theorem (FDT) . Conversely, 
for t ^ tyj they visibly drift and the FDT is violated. As was recently discovered, the drift 
happens in an intermittent fashion [6,7], i.e. through rare, large, and spatially heterogeneous 
re-arrangements, which appear as non-Gaussian tails in the probability density function (PDF) 
of configurational probes such as colloidal particle displacement [8,9] and correlation [10] or 
voltage noise fluctuations in glasses [11]. 

As aging phenomena are similar for a broad class of interactions, we seek a mesoscopic 
description, and assume that intermittent events, for short quakes, are the main source of 
de-correlation in non-equilibrium aging. In the framework of record dynamics [12, 13], quakes 
are irreversible and are triggered by (energy) fluctuations of record magnitude. We show how 
this leads to a description of the configurational autocorrelation function, more specifically, 
the dependence of the shape of its PDF on t, t^, the temperature T and the system size N, 
which resembles observations for colloidal gels [10] spin-glasses and kinetically constrained 
models [15,16]. The model PDF is closely approximated by the Gumbel distributions widely 
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used in the literature [16, 17]. The average and variance are given in close form as a function 
of t/tw, T and N . The average and the PDF, standardized to zero mean and unit variance, 
are in excellent agreement with spin-glass simulations. The agreement is rather poor for the 
variance itself, mainly because pseudo-equilibrium fluctuations are neglected. 

The configuration auto- correlation PDF. - In the model, a set of N binary variables 
defines the system configuration. Without further loss of generality, and with an eye to the 
simulations of the E-A spin-glass model [18], we refer to these variables as spins, and to their 
changes of state as 'flips'. 

Configuration changes are gauged by the number of spins with different orientations at 
times t^ and t^^ + t. This Hamming distance, iJ, is simply related to the autocorrelation C 

by 

C{t^,t) = l-2H{t^,t^ + t)lN. (1) 

Initially, we focus on the probability PH{h,tw + t \ 0,tiu) for H = h at time t^u + t, given 
H = at tyi; , which we write as the average 

oc 

PH{h, i„ + t I 0, t„) - J2 Psis)PHih I s) (2) 

over the conditional probability P/f(/i | s) for H = h given s flips (s = 0, 1, . . . oo). The weight 
function Ps{s) is the probability for exactly s flips during [tuj,tw + t). 

Assuming for simplicity that flips occur at any site with probability 1/A^ leads to the 
master equation: 

PH{h\s + l) = {l-'^)PH{h-l\s) + ^PH{h + l\s) 0<h<N, (3) 
with the initial condition 

PH{h\0)=Sh,o. (4) 

The equation has the formal solution 

PH{h\s)=T'PH{h\0), (5) 

where T is the (bi-diagonal) stochastic matrix implicitly given by Eq. |31 and where the vector 
Ph has elements Ph{0 \ s),Pf/(l | s)...Ph{N \ s). The s dependence of the conditional 
average and variance of Mff(s) and cr|^(s) can be gleaned from the moment generating 
function J2h=o-^H{h \ s)z^, | ^ |< 1- Omitting the details, one flnds 

M/f(^) = f (l-(l-2/A^r) (6) 

and 

N 

^U^) - ^(1 - (1 - 4/iV)^) + —((1 - - (1 - 2lNr). (7) 

Since (Jh{s) <C ^inis) for large N, the r.h.s. of Eq. |2is dominated in this limit by the term 
with index s{h) implicitly given hy h = ^jLh{s). As a consequence, Ph and Ps acquire very 
similar shapes when standardized to zero average and unit variance. 

To calculate Ps{s) we need the probability that i quakes occur between i^, and t^+t and 
the distribution of the number of flips, for short 'size', of each quake. According to refs. [12,13], 
i has a Poisson distribution with average 

niit^,t) = aiN)lii{l+t/t^). (8) 
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(a) Three model PDF's (b) Comparison with Gumbel densities 

Fig. 1 - (a): Three PDF's of the model auto-correlation function C are shown in the main panel, 
for different values of n/. The PDF's are all shifted to zero average. The inserts show the nj 
dependence of the average and the variance, the latter multiplied by A'^. The symbols are from 
numerical evaluations, and the full lines are according to Egs. 1151 and 1161 (b) The model PDF for 
ni = 16.22 is shown together with its best Gumbel approximation, with both curves shifted to zero 
mean and rescaled to unit variance. The left insert shows the linear relation between the model 
parameter nj and the g value for the best Gumbel fit. In the right panel the approximation error, 
(multiplied by 1000), is plotted versus n/. 



The property a{N) cx iV, which removes the N dependence of the exponent A, (see Eg. I14|l . 
arises when the intermittent signal results from independent intermittent processes, stemming 
e.g. from locally thermaUzed clusters [13]. The T independence of a reflects the noise insen- 
sitivity of record dynamics [12], and holds within the low temperature range for which the 
description applies. 

For the quake size, simulations of vortex dynamics [19] yield a near exponential distribu- 
tion. The same form is consistent with the (asymptotically) exponential distribution of the 
energy released [13,20] by intermittent events. Hence, glossing over the integer nature of the 
sizes, we treat them as independent stochastic variables X^, fc = 1, 2 . . . ?, with the PDF 

Px{x)^q{T)eM-q{T)x). (9) 

A temperature dependence of the reciprocal average quake size, q{T), is allowed (but not 
required) by the theory, and is directly observable through the exponent A, see Fig. |21 For 
typographical clarity this dependence is left understood, together with the dependence of ni 
on ln(l -I- t/tw)- 

Considering first the conditional probability for Si flips for a given number i of quakes, we 
note that Si ~ X]fc=i gamma density 

^'s.(x) =<z|^^exp(-ga;). z > 0. (10) 

Averaging the above expression over the Poisson distribution of i, and taking into account 
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that 6(s) exp(— n/) is the probabiUty of no flips (i — 0), one finds 

oo ^ 

Ps{s)=Y,PsM'^exp{-nj)+5{s)cxp{-nj). (11) 

With the variable z = (4gsn/)^/^, this is rewritten as 

Ps{s) — 2njaex.p{—qs — nj)Ii{z)/ z + S{s) exp(— n/), (12) 

where Ii is the modified Bessel function of order one (See e.g. Abramowitz & Stegun, 9.6.10). 

The S{s) term in Eq. ^]will be neglected, since nj is large except for t t^^. With the 
term discarded, the standardized Ps{s) has no T dependence. This is seen, in brief, as follows: 
Using us — nil q and (j| = 2ni/q^ for the average and variance of S, the standardized PDF, 
(JsPs{{s — iJ.s)l'Js), has no q dependence. However, as q carries the model full T dependence, 
the latter disappears as well. Furthermore, due to its similarity with P5, the standardized Ph 
is also independent of T, as confirmed by Fig. |21 

Averaging Eqs.Eland 13 over Ps{s) , reintroducing the time dependence of nj and making 
the approximation — ln(l — 2/N) k, 2/N , one finds 



^JiH{t^.t) = ^{l-{l + t/t^f'^^)), (13) 

where 

The average (macroscopic) form of the correlation function C, is obtained from Eos. 1131 andlTl 
as 

tic{t^,t) = {l + t/t^)^^^\ (15) 

Similar steps lead from Eq. [3 to 

Nal{t^,t) = 1 - (1 + t/t^f^^'^^ (1 - 2A(T) ln(l + t/t^)) . (16) 

Panel (a) of Fig. ^ shows, for three different values of n/ = a(A^)ln(l + t/t«,), the model 
PDF given by Eqs. ^1 [3 and [3 The two inserts show the ni dependence of and Na^ 
from a numerical evaluation of the model equations (circles) and from Eas. 1151 and 1161 with 
the ni{t,tw) dependence reintroduced (fines). 

In standard form (see e.g. ref [16]), the one-parameter family of Gumbel densities is given 
by = r{fy ^MKv - Vo) - e'^y-y°^), with b = ^(^'(ff)) and byo = ln(g) - ^{g), where 

\1/ denotes the digamma function, its derivative and 5 is a real number. Gumbel densities 
empirically describe fluctuations in complex systems [16,17]. Fig.n('6j shows that, except for 
n/ < 1, our model PDF is closely approximated by the Gumbel PDF whose g value minimizes 
the Li distance between the two. The left insert of the figure shows that this optimal g value 
is linearly related to ni as g — O.SOOn/ — 0.185. 

Comparison with simulation data. - The (average) autocorrelation function for spin- 
glasses is well investigated [14,16,21,22]. For t > t^, fJ-c is nearly a function oit/t^, and can 
be fitted by a power-law with a temperature dependent exponent. E.g. Picco et al. [21] found 
an excellent scaling using the variable ln(t^, + i) — In(i^). The autocorrelation PDF for the 
E-A spin glass model nearly follows scaling according to Castillo et al. [15] Chamon et 
al. [16] also consider a kinetically constrained model with trivial statics. In both cases, the 



Paolo Sibani: Mesoscopic statistics and intermittency 



5 




(a) mean of autocorrelation (b) variance of autocorrelation 

Fig. 2 - (a): The mean of the autocorrelation function in an E-A spin-glass of size A*" — 16^ for 
different values of and t is shown by symbols versus ln(l + t/t^), for the indicated temperatures. 
The full lines are power-law fits according to Eq. 1151 The insert shows the temperature dependence 
of the exponent \{T). (b) For the same temperatures, the estimated variances of the autocorrelation 
function multiplied by versus ln(l 4- t/t^) together with la error-bars. The full lines are fitted to a 
shifted and rescaled form of Eg. 1161 as detailed in the text. 



autocorrelation PDF, shifted to zero mean and rescaled to unit variance, is empirically fitted 
to time evolving Gumbcl distributions, which arc numerically equivalent to our model results 
(see Fig. Pi . 

For a more detailed comparison, we simulated the E-A spin-glass [18] on a cubic lattice 
with N = 16^ using an event driven simulation technique [23], whose 'intrinsic' time unit 
corresponds, for large systems, to one Monte Carlo sweep. The data are sampled at 20 time 
points, which are separated by a multiplicative factor of 1.5, with start a,tt = 100 and end at 
t « 2.2 X 10^. Among these points, any ordered pair can be chosen for and + t. For each 
set of physical parameters, 1000 runs are performed with independent noise and couplings 
realizations, producing e.g. 1000 data points for t/tw = 2216, and 20000 points for t/tu, = 1.5. 

For a range of temperatures, the average spin-glass autocorrelation is plotted versus t/t^ 
(symbols). Deviations from t/tw scaling appear for approx. t/t^ < 4 and T > 0.5, as seen 
in the left panel of Fig. |21from the poor data collapse. Away from this parameter region, the 
data are well described by Eq^| (full line). The T dependence of the exponent is shown in 
the insert (circles), where the fit A(r) = -0.25T/Tg is also shown (full line). Tg = 0.97 is 
the critical temperature of the model [24] . The same linear form, and with a similar slope 
coefficient, is found numerically in Kisker et al. [14]. They, however, introduce a small 
dependence of A, which is beyond the present model. 

Summarizing, for low T and not too small t/tw, the model is able to describe the time 
dependences of the average autocorrelation with no free parameters. Furthermore, Eq. 1141 
links A(T) to a linear temperature increase of the average quake size. 

To improve the statistics of the variance and PDF data in the i/im scaling region where a 
comparison with the model is most interesting, we estimate the variance and its error-bar as 
the mean value and standard deviation over the set of variances for all pairs t^^ , t having the 
same ratio i/i^,. Similarly, the empirical frequencies of the C values are calculated based on 
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Fig. 3 - From the left to the right panel, the nj values indicated correspond to t/t^ = 2.3, 7.6 and 
25.6. Within each panel, the scaled and shifted autocorrelation PDF is shown for the model (full 
line) together with four sets of simulation data for temperatures T = 0.15, 0.25, 0.35 and 0.5. 



all data with the same t/tw 

The value of N stands in the model for an (unknown) number of thermally active spins, and 
appears in the autocorrelation variance, which vanishes linearly with This leads to an 

undetermined T dependent scale factor, /i, between model and simulation variance. Secondly, 
and more importantly, the data cannot be fitted without a second, ad hoc, offset parameter 
/2, likely because the de-correlating effect of the pseudo-equilibrium fluctuations is altogether 
neglected. Hence, with respect to the variance, the theory only provides a qualitative de- 
scription, which is captured by the empirical formula Na^^^ {twi't) = fi{,N(T'^[tyj,t) -1-/2)- 
For completeness, the latter is plotted (lines) in the right panel of Fig. |21 together with the 
simulation data with error bars. The parameter /i increases with T within the range 1 — 20, 
and /2 remains close to 1/10, 

By contrast, an excellent agreement between predictions and data for the standardized 
PDF is achieved by a simple adjustement of the vertical scale of the latter. This scale, which is 
undetermined from the outset, is fitted to the properly normalized model PDF. The centering 
and rescaling are done using the data average and standard deviation. The results are plotted 
with la error-bars in Fig.|3 The three panels of the figure correspond, from left to right, to 
t/tw = 2.3, 7.6 and 25.6. In each panel, the data shown are for T = 0.15,0.25,0.35 and 0.5. 
Their collapse confirms the anticipated T independence of the standardized autocorrelation 
PDF. The model predictions (full lines) contain one parameter, a{N), whence it is possible to 
determine one value of n/ by data fitting. This was done for [t/t^ — 25.6) — in the rightmost 
panel — yielding n/ = 78 . Eq. |Slthen gives a — 24, whence ni = 20 and 49 for ~ 2.3 and 
7.6 respectively. 
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Conclusion. - Based on the record dynamics description of intermittency [12, 13] the 
model develops the aging properties of the configuration autocorrelation after a deep quench. 
Its predictions for the average autocorrelation and the standardized PDF are accurate at low 
temperatures and for t > tw Together with allied efforts [19,25-27], the present results 
support the view that record-sized fluctuations are important for aging in metastable glassy 
systems. 
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